Analyzing quantum jumps of one and two atoms strongly coupled to an optical cavity 
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We induce quantum jumps between the hyperfine ground states of one and two Cesium atoms, 
strongly coupled to the mode of a high-finesse optical resonator, and analyze the resulting ran- 
dom telegraph signals. We identify experimental parameters to deduce the atomic spin state non- 
destructively from the stream of photons transmitted through the cavity, achieving a compromise 
between a good signal-to-noise ratio and minimal measurement-induced perturbations. In order to 
extract optimum information about the spin dynamics from the photon count signal, a Bayesian 
update formalism is employed, which yields time-dependent probabilities for the atoms to be in ei- 
ther hyperfine state. We discuss the effect of super-Poissonian photon number distributions caused 
by atomic motion. 



I. INTRODUCTION 

Systems comprised of neutral atoms coupled to a single 
mode of a high-finesse resonator belong to the key model- 
systems in quantum optics In the so called strong 
coupling limit an atom periodically exchanges its excita- 
tion energy with the resonator light field. In this case the 
dynamic evolution is governed by a priori entangled light- 
matter quantum states, namely the combined dressed 
states of the atom-cavity system. Due to the symmetric 
interaction, described by the Janyes-Cummings Hamilto- 
nian Q , information about the state of the system can be 
obtained from two complementary partial measurements: 
In the optical domain, experiments rely on the detection 
of photons emitted from the cavity ; whereas in the 
microwave regime the quantum state of atoms transiting 
the cavity field is detected [l(| • 

Optical cavity-QED systems are attractive for applica- 
tions in quantum information science, e.g. for quantum 
networks. The successful demonstration of, for instance, 
the mapping of the coherent state of a travelin g qu bit 
(a photon) to the atomic state memory qubit as 
well as single photon generation [l2j are recent exam- 
ples of significant progress in controlling the interaction 
of a single atom with the cavity field. For the creation of 
two-particle entangled states, promising proposals rely 
on either applying deterministic protocols II 311 or mea- 
surement induced (probabilistic) projection [l4,ll5j|. 

Not only for the application in quantum information 
science, but also more generally for the investigation of 
strongly interacting atom-cavity systems it is vital to un- 
derstand the spin dynamics of one and especially more 
than one atom simultaneously coupled to the resonator 
field. In the work presented here we concentrate on the 
case of one and two atoms and investigate how maxi- 
mum information about their hyperfine ground state can 
be retrieved from the stream of photons arriving at the 



detector. We outline and detail the identification of opti- 
mal experimental settings such as atom-cavity detuning. 

Random telegraph signals, obtained by continuously 
observing quantum jumps between the spin states of a 
single atom, are analyzed by applying a Bayesian update 
formalism. In complementary experiments (l6l. Il7j. the 
photon number state (Fock state) of a microwave cav- 
ity field is interrogated by a stream of circular Rydberg- 
atoms acting as quantum probes. There Bayesian anal- 
ysis has proven to be a useful method of analysis, 
too. Random telegraph fluctuations are a universal phe- 
nomenon observed in many different fields, including a 
large variety of solid-state systems [lq . 

In order to study two-atom dynamics, we identified 
experimental parameters for which the intra-cavity in- 
tensity depends on the number of atoms in a specific 
spin state. The virtue of the Bayesian method is ev- 
ident in analyzing the corresponding telegraph signals, 
for which the atomic state cannot be unambiguously de- 
duced from the measured transmission signal because of 
technical limitations on the signal-to-noise ratio. 

In our measurements we observe fluctuations in the 
transmission exceeding shot noise, which we attribute 
to thermal motion of the atom. We discuss the im- 
pact of this external dynamics on the performance of the 
Bayesian analysis. 

II. EXPERIMENTAL TECHNIQUES 

A. Setup to trap and transport single atoms 

At the beginning of every experimental sequence, a 
controlled number of cesium (Cs) atoms is transferred 
from a magneto-optical trap (MOT) into a standing wave 
far-off-resonant dipole trap (FORT) with Afort = 1030 
nm and a trap depth of [/fort ~ fc_e x 1 mK. This 
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FIG. 1: Schematic setup of MOT, far-off resonance dipole 
trap (FORT) and cavity mirrors (not to scale). Details on 
the experimental setup and the stabilization of the cavity res- 
onance frequency are given in Q. 



trap acts as an "optical conveyor belt" [19| to trans- 
port atoms into the optical resonator. The fundamental 
TEMoo mode of the cavity has a diameter of = 46 
jim and a length given by the mirror distance of 158 pm, 
the finesse is T = 1.2 x 10 6 . 

A conceptual drawing of the main components is de- 
picted in figure [TJ for details on the cavity-setup and the 
stabilization scheme see ||. To study the atom-cavity 
system, the transmission of a weak probe laser through 
the cavity is detected with a single-photon counting mod- 
ule (SPCM). Using a custom-build time-to-digital con- 
verter, we record - for each photon click - the time since 
the last click, where for our typical count-rates dead time 
effects are negligible. This list of click-intervals is then 
converted into a binned transmission signal by counting 
the detector clicks in each bin time interval Ati,. 

The total detection efficiency for the probe laser light, 
including absorption and scattering by the mirror coat- 
ings, losses at various optical elements along the opti- 
cal path, and the quantum efficiency of the detector, 
amounts to rj — 4.4%, which is a threefold improvement 
compared to our earlier work presented in Q . The main 
challenge was to optimize the separation of probe- and 
stabilization-laser, with typical powers of a few 10~ 15 
and 10" 6 Watts, respectively. In a first step, they are 
separated by their carefully adjusted orthogonal polariza- 
tions. Improved spectral filtering was achieved by replac- 
ing a standard ruled diffraction grating with a volume 
holographic grating, allowing us to omit an additional 
interference filter used before, while still achieving a to- 
tal suppression of the stabilization laser to better than 
lCr 8 . 

The probe laser frequency is set close to the 
\F = 4) — > \F' = 5) transition of the Cs Z?2 hue, where 
F is the total angular momentum quantum number. For 
this transition, the important parameters of the atom- 
cavity system are (3,^,7) = 27r x (13.1,0.4,2.6) MHz, 
where g is the nominal coupling strength for an atom at 
the position of maximum coupling, k is the cavity field 
decay rate, and 7 is the atomic dipole decay rate. Since 
in our setup the birefringent splitting of the cavity res- 
onances is larger than the cavity linewidth, the cavity 
field is always linearly polarized, causing a distribution 
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FIG. 2: Simplified Cs level scheme, (a) An atom in \F = 3) 
is so far detuned from the cavity resonance that it does not 
alter its transmission, (b) If the atom is in \F — 4), it changes 
the transmission, depending on the cavity-atom detuning A ca 
and the coupling strength g. 



of the atomic population over all Zeeman sublevels due to 
photon scattering by the probe laser. Thus the coupling 
strength g given above is obtained from a weighted aver- 
age over all couplings g{mp), based on the steady state 
mf distribution for linearly polarized optical pumping 
[20I l2lT | . With the single-atom cooperativity parameter 
C\ = g 2 /(2K'y) ^> 1, our system is in the strong cou- 
pling regime, where already a single atom significantly 
influences the cavity spectrum. 



B. Nondestructive state detection 



In our system the two long-lived hyperfine ground- 
states \F = 3) and \F = 4) serve as qubit states [22|]. For 
the coupled atom-cavity system we measure this state by 
tuning the cavity close to the \F = 4) — > \F' = 5) transi- 
tion, where only an atom in the \F = 4) state leads to a 
drop in the transmission, while an atom in \F — 3) is so 
far detuned (around 9.2 GHz) that it effectively decou- 
ples from the system and does not influence the cavity 
transmission, see fig. [2] 

The probe laser with angular frequency uj p is initially 
tuned to the resonance frequency of the empty cavity 
lo c = uj p , so when an atom in \F = 4) is inserted into 
the cavity the transmission is reduced to a level which 
depends on the detuning A ca = lo c — co a , where cu a is 
the angular frequency of the atomic \F = 4) — > \F' = 5) 
transition, including the AC-Stark shift induced by the 
FORT potential. To experimentally distinguish between 
an atom in \F = 3) and an atom lost from the trap, which 
both result in the same transmission signal, a repumping 
laser resonant with the \F = 3) — > \F' = 4) transition can 
be applied from the side which brings the atom in \F = 3) 
back to the 1^ = 4) — >■ \F' = 5) cycle. Thus for an empty 
cavity the transmission would remain unchanged, while 
for an atom still present in the cavity the transmission 
would drop again. 
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FIG. 3: (a): Black and grey curves show two single traces of 
quantum jump measurements. The arrows indicate insertion 
and removal of an atom. At the end of the sequence, the 
repumper is switched on again to check that the atom was 
not lost, (b): Ensemble average over 31 single traces. The 
average dwell time is obtained from the exponential fit. 
The averaged transmission level at the end of the sequence, 
when the repumper is switched on, is higher than the initial 
drop, indicating a lower average coupling strength. This could 
be caused by increased thermal motion, a re-distribution over 
different mF-levels, or a combination of both effects. 
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FIG. 4: (a) Normalized one-atom-transmission as a function 
of the cavity-atom detuning A ca . The solid lines are calcu- 
lated for an atom at rest with an effective coupling strength 
of <7eff/(27r) = 8,9, and 10 MHz for the upper, middle, and 
lower curve, respectively, (b) Average dwell time R± 3 as a 
function of detuning. The shaded area is the result of a the- 
oretical model taking motion of the atom into account, and 
the range of values represents our limited knowledge about 
the exact distribution over the Zeeman-sublevels. 



III. SINGLE ATOM SPIN DYNAMICS 

If the state detection technique described above gave 
the same result for an unlimited scries of state measure- 
ments, it would be a perfect projective quantum nonde- 
molition (QND) measurement [23T - I25I ]. assuming the sys- 
tem is otherwise unperturbed. However, in our situation 
the same laser that we use to detect the atomic state (the 
probe laser) can change it via inelastic hyperfine-state- 
changing Raman scattering. An atom in the \F = 4) 
groundstate can thus be transferred to \F = 3) via the 
\F' = 3) and \F' = 4) excited states, thereby undergoing 
a quantum jump [26l - {30j . 



A. Quantum jump rate and transmission level as a 
function of detuning 

To experimentally determine the rate R43 of probe- 
laser-induced transitions from \F = 4) to \F = 3) and to 
identify optimum experimental conditions, we performed 
the following measurement: An atom, optically pumped 
into \F = 4), is transported into the cavity center, caus- 
ing a drop of the cavity transmission, see fig. [3] (a). Since 
no repumper is applied, probe-laser scattering causes a 
spontaneous transition to \F = 3) after some time, visi- 
ble as an instantaneous rise in transmission back to the 
empty cavity level. To check whether the rise in trans- 
mission is really due to a quantum jump and not caused 
by atom loss, the repumper is switched on at the end of 
the sequence as discussed above. 

For each experimental realization, the quantum jump 
occurs at a random point in time, see fig. [3] (a) for two 
example traces. Since the rate of state transitions is time 
independent, the ensemble average plotted in fig. [3] (b) 
reveals the expected exponential curve with the time con- 



stant being the average dwell time i? 43 . 

This average dwell time was measured for a wide range 
of detunings A ca /(2n) = 38. ..410 MHz. For the same 
settings, but with the repumper constantly applied, we 
measured the transmission level T\ , defined as the photon 
count rate with one atom in \F = 4) coupled to the cavity, 
normalized to the empty cavity signal. The results of 
both measurements are presented in fig. [¥] 

In order to describe our measurements with a simpli- 
fied analytical model, we consider a two-level atom at 
rest with the probe-laser being resonant with the empty 
cavity {u p = lo c ). In the weak excitation regime, the one- 
atom-transmission level can be expressed analytically as 

m 

Tx(A ca!5eff )= (7K + ^ )2 + (AcaK)2 . (1) 

The distribution over Zeeman sublevels, thermal motion 
of the atom, AC-Stark shift variations, and other conceiv- 
able perturbations are all accounted for by an effective 
coupling strength g c g. It is defined by Eq. ([1]) in such 
a way that a stationary two-level atom with a coupling 
strength of g e s would yield the experimentally measured 
transmission level. The solid lines in fig. 2] (a) are cal- 
culated according to Eq. (fTJ) with g e g/(2n) = 8,9, and 
10 MHz, and this range of effective couplings describes 
the data reasonably well. We attribute the difference be- 
tween the nominal coupling strength of g/(2n) = 13.1 
MHz and g c g mainly to thermal motion of the atom. 

To describe the measured average dwell times theoret- 
ically, see fig. S] (b), i?43 is calculated as a function of de- 
tuning using the Kramers-Heisenberg- formula [32[. For 
this calculation, one has to treat the distribution over all 
Zeeman sublevels and thermal motion separately, since 
this situation cannot be modeled as a two-level system 
with an effective coupling. The measured data agree sat- 
isfactorily with the theoretical model, confirming that the 
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FIG. 5: (a) Random telegraph signal for one atom coupled to 
the cavity, (b) Bayes analysis yielding po(t), i.e. the proba- 
bility to be in \F = 3). 



best approximation to a projective QND-measurement 
with longest dwell times is close to resonance. A practi- 
cal limitation is that stable coupling was never observed 
for detunings A ca < 2nx 30 MHz^probably due to cavity 
cooling becoming less effective [H, H3] . 



B. Statistical analysis of single-atom random 
telegraph signals 

In the experiments discussed so far, the repumping 
laser was either switched off or its intensity was adjusted 
such that an atom off-resonantly transferred to \F = 3) 
was pumped back to \F = 4} immediately, compared to 
all relevant time scales in our experiment. In contrast, 
for the measurements presented in the following, we de- 
liberately attenuated the continuously applied repump- 
ing laser to a level at which the transfer rate R34 from 
\F = 3) to \F = 4) was comparable to R43. Therefore, 
the resulting quantum jumps occur in both directions on 
a similar timescale of several milliseconds and are thus 
detectable as a random telegraph signal, see fig. [5] (a) for 
an example trace. 

We quantify our knowledge about the atom's hyper- 
fine spin by probabilities assigned to the different atomic 
states. For the following discussion we introduce the pa- 
rameter a to denote the number of atoms in | F = 4) . In 
this section a attains only the values and I, whereas 
the case of two atoms (see section ITV]) also permits the 
value a = 2. Internal state changes of a single atom are 
transitions between the two states a = 0, f , and they 
occur with the rates Rio and i?oi: which are identical 
to R43 and i?34, respectively. Although we imagine the 
transitions to occur randomly and at discrete instances 
of time, the probabilities for the atom to occupy the dif- 
ferent states change in a continuous manner governed by 
the following rate equations: 



dp {t) 
dt 

dt 



-R olPo (t) + R loPl (t) , (2) 

dpo(i) 



The average steady-state probabilities po and pi are ob- 
tained by setting = = and using p + p\ = 
1. The solutions are thus given by the ratio between the 
transition rates: 



Po = 
Pi 



10 



Rio + R01 
Rqi 



R 



10 



-R01 



(4) 
(5) 



The average probabilities and thus the ratio of the rates 
Rio and Rqi can therefore be obtained from photon count 
histograms by the following procedure: Along with the 
telegraph signals, transmission traces for an empty cavity 
(a = 0), and for one continuously coupled atom (a = 1) 
were measured for otherwise identical settings. From 
these three sets of data, normalized photon count his- 
tograms "P(n), V(n\0) and V(n\l) are computed, with n 
being the number of photons detected per binning time 
Ait = 1 ms. Here and for the remaining discussion, V 
always refers to photon-count probabilities, while p a (t) 
indicates spin-state probabilities. 

Since the telegraph signal is expected to represent 
the atomic system jumping between the different states, 
the associated accumulated histogram of photon counts 
should be a weighted fit 



V{n) =poV{n\Q) + {l~po)T(n\l) 



(G) 



of the independently measured histograms V(n\a) for the 
two atomic states. Treating po as a fitting parameter 
yields p — 0.64 and pi = 0.36. 

In order to extract the transition rates, we note that 
the jumping of the atom between two different states 
with different transmission properties causes character- 
istic fluctuations in the number of detection events ob- 
tained in different time bins, n(t) and n(t + t), which 
become visible in the second-order correlation function 
g^(r). Assuming Poissonian count distributions, an 
analysis of the rate equations yields [28| 



9 (2) (r) 



(n{t)n(t + T)} 
(n(t))(n(t + r)} 



<x exp (— (Rio + Rqi)t) for t > 

(7) 



-R w pi(t) + R iPo(t) 



dt 



(3) 



The histogram of the telegraph signal and the func- 
tion are plotted in figure [S] (a) and (b), respectively. 
From an exponential fit of the correlation function, we 
get i?io + i?oi = 50 s _1 , therefore we obtain Ri =40 s _1 
and i?oi = 18 s _1 using p , pi, and Eqs. (|4|) and ([5]). 

In the discussion above, we assumed that the state of 
the atom can be described by the two states a = and 
q = 1 alone, each leading to a Poissonian distribution 
V(n I a) of the photon count rate . For a = this is veri- 
fied by the measurement: The right peak of the measured 
histogram in figure [S] (a) agrees with a Poissonian distri- 
bution of the the same average count rate. Thus the 
state detection for a = 0, with the transmission being 
equal to the empty-cavity case, is essentially shot-noise 
limited and residual frequency or intensity fluctuations 
of the probe-laser can be neglected. 
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FIG. 6: (a) Normalized histogram extracted from 13 tele- 
graph signals of 1000 ms duration each, binned with Atf, = 1 
ms. The solid line is a Poissonian distribution with the same 
maximum count rate, (b) Averaged second order correlation 
function gr' 2 ' for the same set of telegraph signals. The blue 
dashed line is an exponential fit yielding the time constant 
{Rio + Roi)- 1 =20 ms. 



However, comparing the photon count histogram 
V(n\l) with the Poissonian distribution (left peak of the 
histogram in fig. [6] (a)) indicates super- Poissonian fluc- 
tuations. We attribute these mainly to thermal motion 
of the atom: The coupling constant g follows the cavity 
mode function, i.e. g(r) — g^{v) 1 which in turn leads 
to a transmission level ?i(r), depending on the atomic 
position, according to Eq. ([1]). 

In later sections of this paper we shall discuss candi- 
dates for a more complete theoretical analysis of this dy- 
namics. At this point, we pursue a pragmatic approach 
and still extract the atomic transition rates from the cor- 
relation function as stated by Eq. (U), because this rela- 
tion does not rely strongly on the Poissonian character of 
the signal. Furthermore, the exact values of the rates are 
not the main result of this work and do not convey fun- 
damental physical insight, since they are determined by 
the intensities of the probe and repumping laser. They 
rather constitute parameters in the following statistical 
analysis. 

To quantify the knowledge about the atomic state that 
we obtain from the measured telegraph signals, we use a 
Bayesian approach in analyzing the data. The philos- 
ophy behind this approach is that we assign probabil- 
ities to the possible states a = 0, 1 of the atom, and 
acknowledge that these probabilities merely reflect our 
incomplete knowledge about the system, unless one of 
the probabilities is unity. Due to the atomic transitions 
which occur without our direct noticing, the probabili- 
ties of the unobserved system obey the rate equations ([2]) 
and dH) , but since the cavity transmission depends on the 
atomic state, we learn about the atomic state from the 
observed photon count record. 

The probabilities p a are thus calculated step- wise from 
the incremental information obtained in every time bin of 
the measured telegraph signal. Let n(t,) be the number 
of photons detected during the interval [t* — At;,/ 2, ft + 
Atb / 2] , where the binning time At;, is fixed to 1 ms for the 
following analysis. With p a (ti) we refer to the probability 
for an atom to be in the state a in the midpoint of the 
aforementioned interval. Assuming that the atomic state 



probabilities in the previous time binp a (ti_i) are known, 
the probabilities p a (U) are estimated by first evolving 
their values according to the rate equations ((2]) and ((3|). 
In a linear approximation for R x Atb -C 1, where R x — 
max(i?io, Rqx), this leads to 

Po(U) = Po(U-i) + (RioPi(U-x) - RoiPo(ti-i))Atd$) 
Pi(U) = pi(t l _ 1 ) + (i?oiPo(t l -i)-i?i Pi(t l -i))At fl (9) 

where p indicates the unconditional probability. 

Note that the probabilistic description does not imply 
that the atom occupies two different states, but only that 
we do not know which one is actually occupied. This also 
implies that our prediction for the distribution of photon 
numbers n(tj) detected in the i-th time bin has to be 
calculated as a weighted average V(n) = poV(n\0) +(1 — 
pQ)V{n\l). The actually measured photon counts n(ti) 
provide new information, and the state probabilities are 
updated using Bayes' rule of conditional probabilities: 



Pa(U) = p{a\n(U)) 



T(n(U)\a)p a (U) 
T, a Pa(ti)V(n(ti)\a) 



for a — and 1 
(10) 

The conditional probabilities V{n{ti) \a) are extracted 
from the separately measured photon count histograms 
for a = and 1, respectively. Setting the initial prob- 
abilities to Po(Q) — 0, pi(0) — 1, because the atom is 
prepared in 1^ = 4) before being transported into the 
cavity, p a {ti) is then updated stepwise for each time bin. 
In this way, the time-dependent atomic state probabili- 
ties are computed successively for the whole transmission 
trace. 

Figure [5] (a) shows an example trace of a telegraph sig- 
nal to which the Bayesian algorithm was applied. Most 
of the time, the probability po(t), plotted in (b), is close 
to either or 1, while narrow spikes indicate short peri- 
ods of time with less complete knowledge about p a . The 
Bayes analysis provides more definite probabilities than a 
matching of the currently transmitted signal to the state 
dependent transmission rate, because it updates previ- 
ously estimated results and thus accumulates statistical 
significance over time. The optical probing of the sys- 
tem does not, however, prevent atomic transitions from 
taking place, and during such transitions, the Bayes al- 
gorithm faithfully reproduces our inability to determine 
the state of the atom with certainty until a significant 
amount of data has been accumulated which is in agree- 
ment with the new state of the atom. In this context, 
the narrow spikes in the figure illustrate the "willing- 
ness" of the Bayesian update to interpret a few unex- 
pected photon counts as the emerging signal of a change 
of state, while they may be only statistical fluctuations. 
For photon count histograms with negligible overlap, the 
Bayesian algorithm would yield the same result as a sim- 
ple threshold analysis. Its main advantage is that one can 
still extract information about the spin dynamics even 
for a signal where the signal-to-noise ratio prohibits a 
threshold analysis, as will become apparent in the next 
section. 
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IV. SPIN DYNAMICS OF TWO ATOMS 

So far we have presented experiments revealing the in- 
ternal spin dynamics of one atom coupled to the cavity 
mode. Placing two atoms into the resonator leads to an 
effective interaction between them, mediated by the cav- 
ity field |13| , and detecting the number of atoms being in 
a particular state could be used for entanglement gener- 
ation in cavity-QED-systems 



A. Counting the number of atoms in \F = 4) 

In the previous section the atomic state was deter- 
mined from the probe laser transmission. Without 
changing the experimental settings, this is not directly 
possible for two atoms coupled to the resonator. Both 
atoms in \F = 3), i.e. a = 0, will lead to a transmis- 
sion level To = 1 equal to the empty-cavity case. One 
atom in \F = 4) and one in \F = 3) (a = 1) will cause 
the transmission T\ to drop almost to zero, which implies 
that a = 2 is indistinguishable from a = 1. To deduce 
a = 0,1,2 from the corresponding transmission levels 
To, Ti, Tii the experimental settings have to be adapted. 

In the weak excitation limit, two atoms at rest coupled 
with the same strength g to the cavity can be theoreti- 
cally described as a single atom experiencing a coupling 
strength g 2 = V2g- In the dispersive limit (A ca 3> 7), 
Eq. (P) thus yields 



Ti 



1 



T 2 = 



1 



1 



(11) 



1 



for the transmission levels. The level difference 
AT12 = T\ — Ti reaches its maximum value of 33% for 
g 2 /{nA ca ) = l/y/2, where T , Ti, and T 2 are equally 
spaced. In order to examine this theoretical prediction 
experimentally, the transmission level T2 was measured 
alongside the one-atom transmission. Figure [7] shows 
that for two atoms the transmission is lower, but in- 
stead of the theoretically expected value of v2 x 9 MHz 
w 13 MHz, it is compatible with an effective coupling of 
ff2,off ~ 27r x 11 MHz. As a consequence, the measured 
level difference AT12 is at maximum about 20% for a 
detuning of A ca = 2tt x 270 MHz. 

A detuning of A ca in the range of 200 to 300 MHz 
has, however, two disadvantages for studying the spin- 
dynamics of two coupled atoms: Firstly, the difference in 
cavity transmission is quite small compared to the noise, 
and secondly, the average dwell time is close to its 
minimum value for A ca > 2tt x 150 MHz, with a very 
shallow slope towards higher detunings, see fig. 0] (b). 
Closer to resonance, this time is longer, but if two atoms 
are at the cavity center, the transmission levels T\ and 
T2 are almost indistinguishable. 

The level difference AT12 can, however, be controlled 
for a constant detuning A ca by changing <? c ff . This is pos- 
sible by means of our optical conveyor belt, which allows 
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FIG. 7: Normalized transmission T\ (black dots) and T2 (blue 
diamonds) for one and two atoms, respectively. The solid 
lines are calculated according to the effective two-level model 
(TTJ) for one atom at rest with different values for <? c ff, and 
the dashed line shows the theoretically expected two-atom 
transmission for g e s/(2ir) = V2 x 9 = 12.7 MHz. The one- 
atom data is the same as in fig. [4] (a). 



us to transport atoms not only into the cavity center, but 
also to stop the transport at a predetermined distance Ay 
away from it. With g e s(Ay = 0) = 27r x 9 MHz, the cou- 
pling strength as a function of position along the conveyor 
belt axis reads g e g(Ay) = g o ff(0) exp(— Ay 2 /wg). From 
Eq. (frTjl the required distance Ay to achieve AT 12 — 0.33 
is calculated to be 



|Ay(A co )| = w 



\ 



1 ln ( V2g 2 cS (0) 



(12) 



For A ca > 2ir x 280 MHz, AT12 is always at maximum 
for Ay = 0, i.e. at the cavity center. Figure [5] shows the 
calculated level difference AT12 and the quantum jump 
rate R43 as a function of detuning A ca and distance from 
the cavity center Ay. By choosing a lower detuning, the 
scattering rate R43 is reduced, and it is still possible to 
obtain optimal distinction AT12 by positioning the atoms 
away from the cavity center. Empirically we found that a 
detuning of A ca — 2tt x 38 MHz is a lower limit in terms of 
stable transmission traces. The distance of Ay = 21 um, 
at which g e ff/(27r) rs 3.1 MHz, was adjusted for optimum 
distinction of one and two atoms. 



B. Two-atom telegraph signal 

To study two-atom spin dynamics, two atoms loaded 
into the FORT were positioned at Ay = 21 um. At this 
position of around one cavity-waist away from the mode 
center, the coupling strength depends more critical on 
the exact position, therefore those traces were selected for 
which the measured atom-atom spacing was < 2 um. As 
for the one-atom case, the repumper was attenuated to 
a level at which it induced quantum jumps from \F = 3) 
to \F — 4) at a rate comparable with the probe laser 
induced jumps. 
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FIG. 8: (a) Effective coupling as a function of distance Ay 
from the cavity center, (b) Transmission level difference AT12 , 
(c) quantum jump rate R43 as a function of detuning A ca and 
distance Ay. The white solid lines are points of maximum 
AT12 calculated according to Eq. (fT2"l) . 



(a) 




time (ms) 

FIG. 9: (a) Example trace of a random telegraph signal for 
two atoms placed Ay = 21 um away from the cavity center. 
The cavity-atom detuning is A ca — 2ir x 38 MHz. (b) Proba- 
bilities for 0,1, or 2 atoms to be in \F = 4), calculated using 
the Bayes method. 

Figure [§] (a) shows an example single trace of a two- 
atom telegraph signal. For t s» 200 . . . 300 ms, steps cor- 
responding to a = 2 (low transmission), a = 1 (inter- 
mediate level), and a = (empty cavity transmission) 
are discernible, but in general the distinction between 
the levels is not as clear as for the one-atom case. The 
degree of the level separation can be deduced from a his- 
togram extracted from several hundred telegraph signals, 
see fig. [TU] 

This histogram does obviously not show a three-peak 
structure. To quantify the contributions of the transmis- 
sion levels Tq, T±, and T2, we independently measured 
photon count histograms for 0, 1, and 2 atoms coupled 
to the resonator at the same position and for the same de- 
tuning as for the telegraph signals, depicted as solid lines 
in fig. [lOj These were obtained from signals of continu- 
ously coupled atoms, i.e. a sufficiently strong repumper 
was applied. The photon count histogram of the tele- 
graph signal (black line) agrees well with a fit calculated 



1 1 1 1 1 1 1 1 1 




photon count per ms 

FIG. 10: Normalized photon count histogram (bars) of many 
two-atom telegraph signals. The blue, red, and green lines 
are independently measured histograms for 0, 1, and 2 atoms 
coupled continuously to the cavity, respectively. The black 
line is a weighted sum of those three histograms. 



as a weighted sum of the three individual histograms con- 
ditioned on the atomic states. 

The statistical analysis is performed analogous to the 
one atom case, but the set of rate equations now involves 
three atomic states and reads 

^ = -R iPo(t)+R w Pi(t) , (13) 

% = RoiPo(t) - RwPi(t) - Ri 2 Pi{t) + R2iP 2 (tXl4) 
at 

% = Rl2Pl{t) - R 2l P2{t) • (15) 
at 

A transition of an atom from \F = 3) to \F = 4) is only 
induced by the repumper at a rate i? rcp , which is inde- 
pendent of a because the laser is applied from the side of 
the cavity. Thus R12 = i? rep and i?oi = 2i? rop , because 
for the latter case two atoms both in \F — 3) are present. 
In contrast, R21, i.e. the rate that one out of two atoms 
in \F — 4) undergoes a quantum jump to \F = 3), is not 
simply given by 2i?io, because this transition is induced 
by the probe laser, the intensity of which depends on a 
[9| . Theoretically, the jump rate depends linearly on the 
intracavity intensity, thus we expect 

R 21 = 2^R W , (16) 

but this relation is not fixed for the calculation and the 
three rates R21, R21 and i? rep are considered as indepen- 
dent parameters for the calculation. A weighted fit to the 
photon count histogram has two independent fit param- 
eters and yields the steady state populations, which are 
related to the ratio of the three rates. In contrast to the 
single atom case, here it is not possible to the make a rea- 
sonable fit to the correlation function to obtain the sum of 
the rates and thus all three parameters. Instead, we ini- 
tially guess the transition rates and employ the Bayesian 
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update method to extract time dependent atomic state 
probabilities. Then we apply a fit as described below 
to iteratively extract values for the transition rates i?io, 
i?2i and i? rcp , which ensure the optimum agreement of 
the time averaged probabilities with the steady state so- 
lution of the rate equations. 

A good initial guess for the rate Rio can be obtained 
from the transition rate for a single atom placed at the 
same distance Ay away from the cavity center, with no 
rcpumper applied, similar to the measurement presented 
in fig. [3](b). The transmission levels T 2 and T\ are mea- 
sured independently, which yields then an estimate for 
R21 according to Eq. ([16]) , The rate i? rep cannot be 
measured independently, but since the power of the re- 
pumping laser is adjusted such that the transition rates 
from \F = 4) to \F — 3) and vice versa are approximately 
equal, R rcp is set to Rio as a starting value for the cal- 
culation. 

With the initial probabilities po{0) = 0,pi(0) = 
0,^2(0) = 1, the Bayesian algorithm is performed step- 
wise for each time bin as described for the one-atom case, 
yielding probabilities p a (t) . Improved values of the three 
transition rates are now determined by the following it- 
erative, self-consistent method: 

An analytical solution of the rate equations for po{t), 
Pi(t) and P2(t), with the initial conditions given above, 
yields the ensemble-averaged probabilities (p a )(t) with 
the three jump rates as parameters. Averaging over the 
probabilities p a (t) obtained from the analysis of many 
traces provides an experimental result for (p a )(t), which 
can be fitted with the analytical solution, in which the 
rates i?io, R21 and R rep are used as fit parameters. With 
the new values for the rates obtained in this way, the 
Bayes algorithm is applied over again to all experimental 
traces, yielding an updated set of time dependent prob- 
abilities p a (t), which is again averaged to extract the 
rates, etc. The converged set of rates obtained from this 
analysis is 

R 10 = 104s' 1 , R 2 i = 52s- 1 , i? rcp = 45 s^ 1 , (17) 

and the final results for p a (t) for the example trace are 
shown in fig. Fig. |H] (b). The ratio between Rio and 
i?2i obtained from this iterative process does not confirm 
the assumption of Eq. (fTo) , because with T\ « 2T 2 , we 
would expect R\o ~ Ri\- The reason for this discrepancy 
remains unclear at this stage. 



V. DISCUSSION OF STATISTICAL ANALYSIS 

In this section we will address some questions arising in 
connection with the statistical analysis presented in the 
manuscript. First, we will discuss the dependence of the 
Bayesian atomic state analysis on the measurement data 
binning time, which presents interesting questions both 
in the case of Posissonian and non-Poissonian counting 
statistics. Secondly, we will discuss the possible origin of 
the non-Poissonian character of the photon count records 
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FIG. 11: (a) Random telegraph signal with 1 ms binning time, 
(b) Enlarged section of 10 ms showing photon click times. The 
quantum jump occurs at about 24.8 ms. 



and its consequences for our extraction of rate parameters 
and the Bayesian analysis. 



A. Bin size and optimum information extraction 

In the analysis of the one and two atom telegraph sig- 
nals discussed so far, we used binning times of 1 ms. 
Let us recall that the raw-data of the cavity transmission 
consists of a list of time intervals between photon clicks, 
see Fig. [TT] (b) for an example trace. To study some of 
the consequences which a change in time bin size might 
have, we analyzed one and the same set of data using the 
Bayes formalism, but for different bin sizes. 

If long time bins are used, the signal to noise ratio 
in each bin is good, and the count histograms for each 
atomic state become well separated. This implies that 
for long sequences of time, the atomic state probabilities 
will be firmly fixed to values close to zero and unity, 
while the instances where transitions between the states 
occur are not resolved within the duration of a single 
time bin. But this is only true as long as Att < i? -1 , 
with R — max(i?oi ; Rio), because for even longer times 
transitions will occur within a significant fraction of the 
bins causing considerable uncertainty about the actual 
atomic state. 

Going to shorter time bins, the signal-to-noise ratio 
is decreased, and the overlap of the photon count his- 
tograms become larger. Correspondingly, it happens 
more frequently that a less probable, but still possible, 
number of counts in a time bin causes a narrow spike 
in the atomic state probabilities derived from the Bayes 
conditional update rule, where indeed no transition took 
place. This behavior is evident from the spikes in Fig [121 

One would suspect that the additional information 
provided by subdividing data into counts registered in 
the first and second half of every time bin would only 
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FIG. 12: Application of the Bayes algorithm for different bin times. Figures (a) to (c) show the histogram of the telegraph 
signal for a bin time of 1 ms, 100 us, and 10 us, respectively. The transmission signals, generated from the same photon-click 
record, and po(t) are depicted in (d) to (e). 



serve to yield a better estimate of the atomic state, since 
no knowledge is lost by this finer binning of the data. In 
the case where no transitions occur and we aim to detect 
the state initially occupied by the atom in a quantum 
nondemolition manner, the Bayesian analysis indeed be- 
comes independent of data binning size for a Poissonian 
count process. To study this issue in our system, with 
the state-changing rate process occurring simultaneously 
with the probing, we evaluated the one-atom telegraph 
signals for different bin sizes using the Bayesian algo- 
rithm, see fig. [T2J Even for At;, = 10 ps, when there 
is no click in 80% of all bins, the calculated probability 
Po(t) is often close to or 1, although the state probabil- 
ity shows more short spikes compared to Atb = 1 ms. 

To give a single quantitative measure of our uncer- 
tainty about the atomic state, we calculate the entropy 

S=(-J2p^og Pa ), (18) 

a 

where the average (...) is performed over the whole dura- 
tion of all analyzed traces. The entropy is plotted in fig. 
I13l for a range of bin times from 10 ps to 20 ms. The sharp 
rise of S for large bins is due to the high probability in ev- 
ery time bin for an atomic transition to occur. We ascribe 



the increase of S towards shorter bins to the occurrence 
of more spikes in p a (t) , already visible in fig. [12] (d) and 
(e). According to the entropy measure, there seems to be 
an optimum time bin, which is related to the magnitude 
of the quantum jump rates. We recall, however, that the 
entropy (| 18[) is only one of many possible measures of the 
information extracted from the system. If, for example, 
the measurements constitute a component in a feedback 
mechanism, access to data on the shortest time scale may 
yield the better performance with respect to the desired 
goal of the feedback protocol. 

B. Origin and modeling of super-Poissonian count 
distributions 

The existence of an optimum bin time, leading to a 
minimum in the time averaged entropy (|18p . is observed 
both for our experimental histogram data and in simula- 
tions with Poissonian counting statistics associated with 
each atomic state. The case of super-Poissonian counting 
distributions, i.e., distributions with a variance exceed- 
ing the mean value of the number of counts, however, 
presents it own separate problems, and points to more 
elaborate future methods of analysis. 
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FIG. 13: Time- and ensemble averaged entropy S as a func- 
tion of binning time At;,. 



Wc already commented on the apparent extra fluctua- 
tions in the light transmission signal being possibly cor- 
related with the atomic motion between sites exhibiting 
different coupling strengths to the cavity mode, corre- 
sponding to different transmission levels. This suggests 
an extended model, where the state with no atoms cou- 
pled to the field (a = 0) is retained as a single state, 
while states with a =1,2 are split according to an ex- 
tra position label, attaining a number of different values. 
If for example a single atom can reside in two locations 
leading to two different Poissonian transmission signals, 
the long time averaged photon count distribution will 
be a weighted sum of these distributions, while the count 
number correlation function within an experimental trace 
may reveal the transition rates between the atomic loca- 
tions, equivalent to our analysis of internal state transi- 
tions in sec. IIII Bl This is an appealing and very likely 
explanation of the broadened histograms, and it points 
to an interesting problem for our previous analysis. 

If the super-Poissonian fluctuations in our counting 
histograms are caused by atomic motion between states 
with different Poissonian signals, counts in close lying 
time bins, where atoms have not yet moved, should be 
correlated. This implies that the Bayesian update is no 
longer a Markov process, where the updated probabili- 
ties depend only on the most recent value and the lat- 
est measurement result, but also knowledge of previous 
counts should be applied to extract maximum informa- 
tion about the atomic state. This effect may have signif- 
icant consequences for very short time bins, where each 
bin offers a low signal-to-noise ratio, but where correla- 
tions between bins may be strong. We have analyzed 
our experimental records, and we indeed find such cor- 
relations, but because of limited statistics these findings 
could not be incorporated quantitatively into our anal- 
ysis. This does not imply that our previous use of the 
Bayes update formalism produces erroneous results, but 
it should be noted that it represents the update based on 



a restricted access to (or memory of) the measurement 
data, and hence it provides a non-optimal estimate of the 
atomic state based on incomplete information. 



VI. OUTLOOK 

We have shown that a Bayesian analysis of experimen- 
tal transmission signals from a cavity containing one or 
two atoms provides a high degree of certainty about the 
atomic state. We have demonstrated how rates of the 
atomic processes can be fitted to the data, and we have 
discussed possible physical explanations of noise in the 
data beyond the predictions of simple models. 

A natural next step would be to use a more complete 
model, including the larger number of position states and 
internal states of the atoms. We recall that already for 
the simplest model with only internal state dynamics, 
finding the parameters is not a trivial task, but ad hoc it- 
erative procedures have allowed the identification of con- 
sistent sets of parameters used in our present analysis in 
this paper. 

It will put stringent demands on the reproducibility of 
large data sets to make a reliable fit to more advanced 
models, but we wish to conclude this paper with a brief 
mentioning of a promising systematic theoretical data 
analysis that can be applied to such data in a future 
more elaborate treatment: the Hidden Markov Model 
(HMM) (3f| . We indicated that there is a possible physi- 
cal mechanism responsible for the fluctuations and for the 
temporal correlations between count signals. In this way 
we point at an underlying Markovian model, where the 
atoms perform transitions between different internal and 
position states, and for each of these states, the coherent 
light field is transmitted with a definite transmission co- 
efficient, and counting statistics are Poissonian with no 
temporal correlations. This is, indeed, a physical realiza- 
tion of a Hidden Markov Process in statistical modelling 
of time series, with applications in insurance, finance, 
speech recognition, image analysis and many other fields, 
where a single series of data is mathematically modeled as 
the outcome of a system undergoing transitions between 
(hidden) states, each giving different data characteristics. 
In their most advanced forms, HMMs only assume the 
transitions between the hidden states to be Markovian, 
i.e., the state populations follow a transfer matrix of dis- 
crete or continuous population changes, while the signal 
can have any state dependent probability distribution. 

Our problem belongs to a narrower class with continu- 
ous rate equations (with unknown rates), and it is plau- 
sible to assume Poissonian count statistics parameterized 
by a single parameter for each atomic state. This case is 
treated, e.g., in [36|, and the problem of estimating the 
transition rates among a family of N states and the N 
photon transmission rates from the data is solved by an 
iterative variational application of the maximum likeli- 
hood principle. In a genuine HMM, the number of states 
N is not known, and one merely attempts to fit the data 
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with different candidate numbers of states. For an appli- 
cation to our problem, we are guided by the physics, and 
after a successful fit, we would request that the states 
identified should have the properties corresponding to a 
few position states for each of the internal state a = 1, 2 
cases. I.e., they should occur in groups with similar pho- 
ton scattering rates, and certain transition rates should 
be very small or vanish. 

In addition to an extended model for the analysis, we 
aim at improving the experimental conditions, such that 
the super-Poissonian noise is less pronounced. Since we 
attribute these fluctuations mainly to thermal motion of 
the atom, increasing the stability of the coupling strength 
requires a tighter confinement of the atom. This could be 
achieved by employing cavity-mediated cooling forces (3ol 
|34| , Raman cooling 4] or additional trapping potentials. 

The rate at which information about the atom-cavity 
system can be acquired is ultimately limited by the pho- 
ton flux arriving at the detector. The most important 
obstacles for further enhancement of the detection effi- 
ciency are losses from the cavity-mirror coatings and the 



limited quantum efficiency of the SPCM. Employing ho- 
modyne or heterodyne detection would permit the use of 
detectors with a quantum efficiency close to 100%. The 
former problem could be solved by using a more open 
cavity configuration, where the transmission coefficient 
is significantly larger than the losses. 

Advancements both in terms of experimental condi- 
tions and statistical analysis could finally lead to the 
development and implementation of quantum feedback 
techniques for the preparation, stabilization and error 
correction of non-classical quantum states [U [H| . 
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